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SUMMARY 


In this study, calculation of longitudinal and lateral- 
directional aerodynamic characteristics of airplanes by the VORSTAB 
code is examined. The numerical predictions are based on the 
potential flow theory with the corrections of high angle-of-attack 
phenomena; namely, vortex flow and boundary layer separation 
effects. To account for the vortex flow effect, vortex lift, vortex 
action point, augmented vortex lift and vortex breakdown effect 
through the method of suction analogy are included. The effect of 
boundary layer separation is obtained by matching the nonlinear 
section data with the three-dimensional lift characteristics 
iteratively. 

Through correlation with results for nine fighter 
configurations, it is concluded that reasonably accurate prediction 
of longitudinal and static lateral-directional aerodynamics can be 
obtained with the VORSTAB code up to an angle of attack at which 
wake interference and forebody vortex effect are not important. 
Possible reasons for discrepancy at higher angles of attack are 
discussed. 
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1 . INTRODUCTION 


Currently there is a strong interest in the technical community 
in longitudinal and lateral-directional aerodynamics at high angles 
of attack for fighter aircraft* Maneuvering flight at high angles 
of attack may induce boundary-layer separation and vortex-separated 
flow over aerodynamic surfaces. These flow phenomena affect not 
only the individual aerodynamic surfaces on which these phenomena 
occur but also those in close proximity. As a result of an adverse 
pressure gradient, the well organized vortex flow may break down, 
significantly changing the aerodynamics of the vehicle. These flow 
patterns are now well known qualitatively for fighter configurations 
from extensive testing in wind and water tunnels. However, 
development of quantitative prediction methods has not kept pace 
with tunnel testing. 

Eventually, a theoretical method for predicting aerodynamics at 
high angles of attack would be ideally based on Navier-Stokes 
solutions with appropriate turbulence modeling. However, for 
applications to preliminary design, this is not feasible at the 
present time due to inadequate computer resources. Therefore, a 
practical approach would be one involving panel or panel-like 
methods with corrections for high angle-of-attack flow features. 

For example, one currently available computer code for complete 
aircraft configurations is the PAN AIR (refs. 1 and 2). This code 
is based on the solution of the Prandtl-Glauert equation and may 
include boundary layer corrections. However, its applications to 



configurations with extensive boundary-layer separation or vortex 
flow have not been demonstrated. It should be noted that 
symmetrical vortex flow without breakdown on simple configurations 
has been dealt with successfully in reference 3 based on a panel 
method. Similarly, the VSAERO code (refs. 4 and 5) was developed 
mainly for symmetrical flight conditions with vortex flow and 
boundary layer correction at moderate angles of attack. A 
comprehensive method for lateral-directional aerodynamics with 
aer ©elastic effect is that of the FLEXSTAB (ref. 6). However, it 
was valid only in the region of linear aerodynamics, and hence at 
low angles of attack. Recently, the VSAERO code has been extended 
to calculation of stability and control characteristics of airplanes 
(ref* 7). Again, applications to high angle-of-attack aerodynamics 
have not been reported. 

To remedy the inadequacy of existing methods, in particular in 
the area of nonlinear lateral-directional aerodynamics, the VORSTAB 
code was developed for vortex-dominated configurations (ref. 8). 

The effect of vortex lift is included through the method of suction 
analogy. To account for the effect of vortex breakdown, 
experimental data were used through a theoretical correlation 
parameter. The code was later extended to treat conventional 
fighter configurations by using nonlinear section data for the 
effect of boundary-layer separation (ref. 9). The nonlinear section 
data could be experimental or theoretical. This was done to avoid 
extreme difficulty in theoretical prediction of three-dimensional 
viscous separated flows. 



In the present investigation, extensive applications of the 
extended VORSTAB code to fighter configurations in nonlinear 
longitudinal and lateral-directional aerodynamics will he 
reported* Both theoretical and empirical methods used in the code 
will be described* Numerical results for several current fighter 
aircraft will be presented* 

2. THEORETICAL APPROACHES 
2.1 General Concept 

The main flow features accounted for in the VORSTAB code are 
illustrated in Figure 1. The flow model in the code is the 
potential subsonic flow, with corrections for the vortex flow effect 
and the boundary layer separation* The potential flow calculation 
is based on the small disturbance approximation of the gas dynamic 
equations which result in the Prandtl-Glauert equation* There are 
two kinds of flow singularities used in the program: horseshoe 

vortices for lifting surfaces and vortex multiplets for the 
fuselage* The strength of the singularities can be obtained by 
satisfying the boundary conditions on the fuselage and lifting 
surfaces* To properly account for the leading-edge singularity of 
pressure loading in the linear theory, and hence the leading-edge 
thrust, the Quasi-Vortex-Lattice Method (QVLM) is used (ref. 10). 

The lifting pressure distribution can also be correctly calculated 
by this method. 
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In addition to potential flow, the vortex flow also 
significantly affects the high angle-of -attack aerodynamics. The 
characteristics associated with vortex flow include vortex lift, 
action points, breakdown effect, etc. The suction analogy is used 
to compute the leading- and side-edge suction forces generated by 
the vortex flow. The action point is predicted by applying the 
linear momentum principle of fluid mechanics. Vortex-breakdown 
angles of attack and progression rates are obtained from semi- 
empirical formulas derived from analysis of experimental data. 

Another phenomenon limiting the high angle-of-attack flight 
envelopes for high performance aircraft is the boundary-layer 
separation. To account for this effect, sectional nonlinear data 
are used iteratively to account for the viscous effect on the 
lifting surfaces. In the following, methodologies used in the code 
will be summarized. 


2.2 Potential Flow Theory (Quasi-Vortex-Lattice Method) 


As indicated earlier, the present potential flow method is 
based on the solution of the Prandtl-Glauert equation: 

„ . M 2) £i + £t + ii* - o 

* 3x 2 3y 2 3z 2 


( 1 ) 


where $ is the perturbation velocity potential. The solution is 
represented by vortex distributions. The thin wing approximation is 
used throughout. 
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2.2.1 Boundary Condition 


The boundary condition for the Prandtl— Glauert equation on the 
wing is that the normal velocity component to the wing surface 
should be zero. Assume that the wing surface can be described as 


z = z (x,y) 


( 2 ) 


Therefore, a unit normal vector on the wing surface can be defined 


as 


n 


3z 3z 

c t -— £• | + £ 

3y J 


3x 


3Z €\ 3Z 


3x 

The total velocity vector is 

^ = (V cosa + u)i + vj + (w + V sina)k • 
00 00 

Applying the tangency condition 
^ • n = 0 


(3) 


(4) 


the boundary condition becomes 


dz 

< 

3x 


3z 

< 

3y 


-V cosa -r- c - v - tt - + V sina + w = 0 


(5) 


where u3z /3x has been ignored as a second-order term. For a wing 
c 

with dihedral (<f>) and twist (a tw ), it can be shown that (ref. 11) 


dz 


„ -sina + -j — cosa 

3z tw dx 0 tw 

c _ _ ^ 

3x 


dz 

cos<t>(cosa tw + ^ sina tw ) 


(6a) 
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dz 

Q 

sitxd) cosa + -= — cosd) 

tw dy 2 r 

dz c 

cos<f> cosa - “ — sin<b 

tw dy 9 r 


(6b) 


where the (X 2 , y 2 ) coordinate system is illustrated in Figure 2* 
Basically, it is a coordinate system based on the local wing chord 
plane in a nonplanar configuration* For a cambered wing with small 
spanwise slopes (dz c /dy2)> it can be shown from Equation (6b) that 

3z c 

Sy- - tan $ < 7 '> 

Therefore, the wing boundary condition becomes 

w v 

V - •— cos(j) - *= 7 — sin <b 
n V Y V Y 

OO 00 

9z c 

= coscj) — — cosa ™ cos(j) sina (8) 

With the presence of a fuselage, the induced velocity on the 
wing consists of two terms: the effect of the wing and the effect 

of the fuselage. The boundary conditions on the wing and fuselage 
become, for the longitudinal case, 

Wing 

9z c 

V + V = cos d> — cosa ■*- cos <b sina (9) 

nw nf Y 9x Y ' 


Fuselage 


3 *f . M w dR . 

-f — = cosa -j— - sina cosQ 
3r 3r dx 


( 10 ) 


where R(x) is the fuselage radius. The fuselage coordinate system 
(x, 9, r) is illustrated in Figure 3. 
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For an airplane having angular motions with rates p, q, and r, 
the induced velocity at (x, y, z) due to these angular rates is 

-(■-tp + jfq - £r) x (lx + jy + itz) 

= £(-qz - ry ) + £(rx - pz) + &(py + qx) (11) 

Adding a sideslip velocity, -f3V , to the total velocity, the latter 
becomes 


^ = (V^cosa - qz - ry + u)i + 

+ (v + rx - pz - BV^)* + &(w + py + qx + Vasina) (12) 

Substituting Equation (12) into Equation (4) and using Equation (7), 
the boundary condition for the angular motion becomes, after 
removing steady symmetrical condition and second-order terms, 

Wing 

V = V + V _ = cos <j> - — sin* 
n nw nf V V 

00 00 

= -Bsin* + sin* - cos* 

OO 00 

= -Bsin* + (-||“)(-|2-)sin* - ('ll - ) \ (.7 cos* + z sin*) 

00 OO 



(13) 


OO c 

In the current version of the code, the effect of pitching (q- 
term) has not been implemented. On a circular fuselage, the 
condition is 


3$ 


f 


3r 



3sin6 - (|b_) (lysine 
00 


(14) 
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To satisfy these boundary conditions, two types of flow 
singularities are used in the program: horseshoe vortices 

representing lifting surfaces and vortex multiplets for the 
fuselage. To sum the effect of horseshoe vortices, the QVLM 
methodology is used to account for mathematical singularities of the 
square-root type at the leading edges and the Cauchy type. A two- 
dimensional formulation will be used to illustrate these points in 
the next section. 

For the fuselage effect, G. N. Ward's vortex multiplets are 
distributed along the fuselage axis (ref. 12), The velocity 
potential is given by 


$ f (x, 0,r) - 


4tt 


I { 


n 


cos n0 
sin n9 


} J ft tx - g + /(x - g) 2 + 8 2 r 2 ] n 

x p o n 7 ( 77 - ~~1 o 2 2 

it r /(x - g) + 3 r 


f„ce>d 5 


(15) 


where cos n0 is for the longitudinal case and sin n0 is for the 
lateral-directional case. When Equation (15) is substituted into 
Equation (10) or (14), it appears that the latter can be decomposed 
into Fourier components if the wing effect, is also decomposed 

likewise. Let 



r rCOS n0 

4, a nisin n0 
n=l 


Then, the Fourier coefficients are given by 


a 

n 


ix 9 $ 


_ £ r e W) r cos n9 

tt l ^9r ^sin n9 


}d0 




(16) 


(17) 
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This concept is implemented by calculating the normal velocity on 
the fuselage circumference induced by the wing vortex distribution 
at a number of fuselage stations and performing numerical 
integration of Equation (17), This method of Fourier component 
matching in satisfying the fuselage boundary condition has the 
advantage of requiring a much lower number of flow singularities in 
the solution process. 

It should be noted that when a combination of surface and axial 
vortex distributions is used, the wing loading may become too low at 
the wing-fuselage juncture due to the horseshoe-vortex 
discretization of wing vortex distributions. This problem is solved 
in the present code, under the symmetrical loading condition, by 
eliminating the inboard trailing vortices associated with the root 
vortex strip. This idea is still correct even if the fuselage is 
absent. If this is not done, the root vortex strip will exhibit a 
loading similar to that near the wing tip. However, this step is 
not needed for the antisymmetrical case in the lateral-directional 
motion. 

Detailed expressions for the induced velocity vector to be used 
in the boundary conditions can be found in reference 10 for lifting 
surfaces and reference 13 for a body. Note that the wake of the 
lifting surfaces is assumed to be flat. 
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2. 2. 2 Calculation of Leading-Edge Suction 


To illustrate the essential idea of QVLM, consider a thin 
airfoil in a two-dimensional flow. The induced velocity at x on the 
chord line is 


-C.. 0 , -_4_r ^ 


(18) 


The integrand in Equation (18) possesses a Cauchy singularity at C = 
x. In addition, y(£) has a square-root singularity at £ = 0. To 


eliminate the square-root singularity, the ^coordinate is 
transformed to a 9-coordinate through the following relation: 

(1 - COS0) 
x 2 

f _ (1 - COS0' ) 

2 “ 


(19) 

( 20 ) 


Therefore, Equation (18) becomes 

/ fl x = J. / y( 9 1 )sin9 *d8 * 
2 tt ^ cos 9 - cos 9 ' 

Let 


( 21 ) 


g( 9) = y(9)sin(0) (22) 

Since y(0) has a square-root singularity at 9=0 and sin 9 vanishes 
at 9 = 0 as the square root of x, the square-root singularity 
of y(0) is eliminated by the factor sinQ. Equation (21) can be 
written as 


w(0) 



Equation (23) can now be reduced to a finite sum through the 
midpoint trapezoidal rule as follows: 


( 23 ) 



i = 0 


(24) 


w(x i ) = W 


Vk (1 ~ 

x i " \ 


1/2 


NC 


i 0 


x, at 0. = iir/N 
x 1 

x k at 0 k = (2k - 1 )tt/2N 

where C is related to a f requently~used leading-edge singularity 
parameter C f 


C T = lim u(x)/x = lim y y(x)/x = ^ (25) 

x-K) x*0 

and can be computed once are obtained. 

In a three-dimensional flow, the wing surface is divided into 
vortex strips. Over each vortex strip, the vortex integral is 
reduced to a finite sum in exactly the same manner as in the thin 
airfoil case described above with the leading-edge singularity 
parameter C calculated from the following relation (ref. 10): 


NC/tan 2 A^ + f$ 2 = l 


induced upwash at £.e. 


Bz 

-[cos<f> — cos a - cos<f>sina] 

The sectional leading-edge thrust coefficient 


(26) 

(c t ) is then given by 


2 2 2 

ttC /I - M cos A. 
_ «> Z_ 

C t 2cosA^ 


(27) 


2.2.3 Pressure Distribution and Total Force and Moment Calculation 

In the linearized theory of both compressible and 
incompressible flow, the wing pressure coefficient .(C ) is related 
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to u in the expression 


C « -2 ~ 
p V 

^ rw 


( 28 ) 


In the thin wing theory, the vortex density function Yy 
equals -2 Therefore, the lifting pressure is given by 

oo 

AC = 2y 

P y 

At any angle of attack, the pressure difference can be shown as 
(ref. 11) 


(29) 


AC’ = 
P 


cos(a + a - 6 )• 2y 
tw c y 


(30) 


dz 

“1 c 

where 6 , equal to tan (-- — ), is the angle of the camber slope of 

c dx 2 

the surface relative to the chord line* 

For a cambered wing, additional lifting pressure will be 
generated from the interaction of a freest ream component with the 
streamwise vortex density y x . Adding this component of lifting 
pressure to AC^, the lifting pressure coefficient becomes (ref. 11) 


AC TI = cos (a + a - 6 ) • 2y - 2y sina sin($ + <f> ) 
p tw c y x y 

, dz 

where <p is the dihedral angle, ^ = tan ^ 


(31) 


With the pressure distribution calculated, the sectional force 
coefficient in potential flow can be calculated by chordwise 
integration as 


+ —Pi 1 ♦ & 2 + (^ 2 ] 1/2 ^ < 32 > 

X le 
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« 4 „ ■ 7 / Xte ^[- + — ]/i‘ * c£) 2 + &N n ** (33) 

X £e 

In an attached flow, c A must be reduced with the leading-edge 

a, p 

thrust. The total force coefficients are calculated by spanwise 

integration of sectional force coefficients as 

b/2 


%P 


Xp S 


/ 

0 


C A,p Cdy 


= - / 
5 J 

0 


b/2 


c , cdy 
d,p 


(34) 

(35) 


The sectional pitching moment coefficient is 


m,p 


2 

c x 


, X. 3 z oz 0 oz - , 

l -l ^)/l» + (^) 2 + (-sr) ] ' o* 


3z 

< 

3x 


3z 


•3y 


(36) 


He 


where Ax and Az are the moment arms of the horseshoe vortex 
elements. Therefore, the total pitching moment coefficient in 
potential flow for a lifting surface is 


m,p 


~I 

Sc 0 


b/2 


m,p 


c dy 


(37) 


The pressure distribution on a body of revolution is computed 
by the following second-order formula: 

C ,_v = -2u - [(1 - M 2 )u 2 + v 2 + w 2 ] (38) 

p(f) 00 

where u, v, and w are nondimensional perturbed velocity components 
referred to the freestream velocity. Converting the velocity 
components u, v, and w to cylindrical coordinates. Equation (38) 
becomes 

G ^ ^ = 1 - (1 + u) 2 + M 2 u 2 + (cosa ~|> 2 + (“ Yfi ~ sina sin0) 2 

(39) 
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where the dR/dx-term is obtained with Equation (10) for the normal 
induced velocity on the body surface. The fuselage sectional normal 
force coefficients can be calculated by integration along the 
circumference as 


j 2 tt tt 

c n(f) = " r l C p(£) cos9 rd6 = ~ 2 f C p(f) cos9 d9 
By integrating the sectional normal force coefficients 
fuselage, the total normal force can be calculated: 


(40) 


along the 


N(f ) 


= i ! 


S ' rc n(f) dx 
fZ 


(41) 


x c 


where x^ is the fuselage station behind which the potential flow 
does not exist (ref. 14). 

Similarly, the fuselage moment coefficient can be obtained from 


"m(f ) 


i x 

= -^/ F 


Sc x. 


rc ,-*xdx 
n(f ) 


(42) 


f Z 


2.3 Vortex Flow Theory 

2.3.1 A Generalized Suction Analogy 


Experiments show that for a low aspect-ratio wing with sharp 
leading edges, the flow always separates from the leading edge and 
rolls up into spiral vortex sheets and reattaches inboard of the 
vortex sheets. A typical flow field of leading edge vortex is shown 
in Figure 4. Due to this kind of flow separation, the pressure 
distribution on a low aspect— ratio wing is different from that given 
by the thin wing theory. In the method of suction analogy (ref. 15), 
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if the flow over the vortex sheets reattaches, the vortex lift on 
the wing is assumed to be the same as the leading edge suction force 
in attached flow. Such vortex separation also occurs along the tips 
(ref. 16). Therefore, it is essential to be able to predict 
accurately the leading-edge and side-edge suction coefficients. 

To calculate the sectional leading-edge suction coefficient 
(c note that the suction force is the force which is normal to 
the leading edge, while the leading-edge thrust is the force in the 
x-z plane. The leading-edge thrust coefficient is calculated with 
Equation (27). For a plane wing, the sectional suction coefficient 
(c g ) is related to the thrust coefficient (c t ) as 

c g ■** c t /cosA£ (43) 

For a cambered wing with dihedral (<}>), the relation is more 
complicated. The latter was derived in Reference 11 and is given as 


C s C t 9z dz 


9z 0 dz 0 ,->■ 

* * (TT? + krl! If l 


where 


and 


{[■ 


z 


9z 9z 


9x dy 




2 , 1/2 


dz o 3z r 

tanA 0 + -»-= 
dy 9x Z 9y 


(44) 


(45) 


dz 


JU 2 


|ii| - /I + tan\ + (^i) 


(46) 


The side-edge suction coefficient will be given later, as it is 
an essential parameter in lateral-directional aerodynamics. 
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In addition, when the leading-edge vortex system passes over an 
area downstream of the leading— edge tip, additional vortex lift is 
developed. This additional lift is called the augmented vortex lift 
(ref. 16). To derive an expression for the latter, consider the 
geometry in Figure 5. If the total vortex lift produced along the 
strake leading edge is equated to the force due to momentum 
transfer, then 

-/ J pV^c s c<ly » /p v v (^ v ' 

- -/ 2p VVav dt) 


It follows that 


av 


4V A 
v 


IT 1 c s cdy 


(47) 


where Z is the leading-edge length over which the vortex lift is 
generated. If the augmented vortex lift (F a ) is assumed to be the 
force produced by momentum transfer over a length (c) as shown in 
Figure 5, then 

F a ' C L(aug) 7 P V “ S “ 2 ’V V v“av> 

' 7 ! / V dy 

Therefore, the augmented vortex lift coefficient is given by 


= £ C T 


J L(aug) Z L(v,£e) 


(48) 


where 


°L(v,jte) - i/ c s cd y 


(49) 
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The sectional aerodynamic characteristics due to the leading- 


edge vortex are given as follows: 


3z 


3z 


3z 


~ O 0 1/0 

= sinot + cosa)/[l + + Gj^-) 2 ] 17 


'l, vie s ^ 3x 


3y 


(50) 


C d,vJle C s v 3x 


3z 3z 0 3z oi/o 

(" — cosa + sina)/[l + (—•) + (-g~) ] / (51) 


C 3z 3z ry 3z ry 1 try 

'* 1 * )/[l + + (-3—) ] 7 < 52 > 


c „ = - — f.Ax n + Az rt ^ 

ra,v£e c v £e ie 3x 


Therefore, the total aerodynamic characteristics due to the leading- 
edge vortex are 



0 

b/2 


L , v£e 

= s' 

l C *,v*e Cdy 

(53) 


0 

b/2 


D, vJte 

6 . 

= s 

/ c, „ cdy 

0 d,vle 

(54) 


2 

r b/2 2 

(55) 

m, le 

Sc 

l C m,vle C dy 


2.3,2 Effect of Leading-Edge Radius on Vortex Separation 

Kulfan assumed that on a slender wing the leading-edge vortex 

separation starts at an angle of attack (a g ) at which the leading- 

edge drag equals the leading-edge thrust (refs. 17 and 18). This 

condition has been shown to be consistent with the angle of attack 

at which the leading-edge laminar separation first occurs (ref. 19). 

To calculate a , consider a cambered wing. The sectional leading- 
s 

edge suction coefficient can be written as 

c = K(sina + a )^ (56) 

s o 
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where K is a function of geometry and Mach number, and a is the 

o 

sectional angle of zero lift. For a wing, the sectional leading- 
edge thrust coefficient was given in Equation 27. The leading-edge 
singularity parameter at any a can be written as 

C = K’ (since + a Q ) (57) 

Therefore, the leading-edge singularity parameter at a Q is 
obtained as 


C, = C(sina + a )/(sina + a ) (58) 

1 s o o 

The starting angle of attack (a g ) of the leading-edge vortex 
separation can be obtained by equating the leading-edge drag to the 
leading-edge thrust. Using the expression for the leading-edge drag 
from reference 20, it is obtained that 


cosA 


✓ . \2 ,,2 2 . x 1/2 

„ (sina + a ) (1 - M cos A n ) 

* JL r 2 S ° °° & _ (zq) 

2 2. x 1/2 2 ^ , . . v 2 cos A, V ' 


(1 - M cos A n ) 
00 £ 


(sina + a ) 
o 


where r Q is the leading-edge radius. It follows that 

, sina + a r , 9 ? i /? 

a “ sin [± — - — 1 -- (2 — ) ~cosA /(l - M cos A ) - a ] (60) 

s C c Z 00 X o 

With a g calculated, the sectional thrust coefficient at a > a g is 

given by 


c = (ir/2)C?(l - M^cos 2 A ) 1//2 /cosA 
t 2 00 Z x. 


(61) 


where 


C„ = C[sin(a - a ) + a ]/(sina + a ) (62) 

2 so o 

This sectional thrust coefficient is converted into the vortex lift 

coefficient through the suction analogy. From Equation (62), it can 
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be seen that the effect of a rounded leading edge is to decrease the 

vortex lift at an a > a . 

s 

2.3,3 Concept of Vortex Action Point 


In developing the concept of the vortex action point, the 

linear momentum principle of fluid mechanics is used. It is assumed 

that the velocity distribution around the vortex in the longitudinal 

plane parallel to the freestream is similar to that in the cross- 

flow plane (Fig. f )• If a control surface o is taken through the 

vortex center as shown in Figure 7, the vortex force should be equal 

to the vertical component of the force due to the momentum transfer 

through the control surface a. As shown previously, the vortex 

1 2 

force acting on the wing is pV^ cc g • It follows that the force 

acting on the control volume is 


"I pV o 


c c 


= / pV (V • dx) = / (pV 


pV„ 


)V 


dr 


in 


out 


To find an average V, let 

|V| = v z = fv w 

where f is a constant. Therefore for a unit span 


(63) 


(64) 


/ (pV - pV )V • dr = -fV [/ pV dt + / pV dr] 

Z in Z out Z in Z out 

= -fV 2p(fV )r (65) 

oo oo v 

where r y is the location of the vortex action point from the leading 
edge. From Equations (63) and (65) it is found that 
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(66a) 


1 

r a 

V 4f 

In the program, f is determined to be 0.5 for subsonic flow (ref. 
11). Therefore, the action point location (r y ) is 

r v = c c g (66b) 


2.3.4 Effect of Vortex Breakdown 


The effect of vortex breakdown is one of the important factors 
affecting the aerodynamics of high performance aircraft. Since 
there is no comprehensive theory to predict the breakdown location 
and the residual vortex strength after breakdown, a semi— empirical 
formula derived from a least-square analysis of available data is 
used in the program. 

Lamar (ref. 16) observed that for delta wings, the angles of 
attack for vortex breakdown at the trailing edge (ot^) is related to 
the leading-edge suction distribution (c g ). Let 


- - 2 

c = c c /c sin a (67) 

s s 

where c is the mean aerodynamic chord (MAC). Since c is 

s 

2 *- 

proportional to sin a, it follows that c g is a function of planform 
and Mach number* From a least-square analysis of Wentz’s data (ref* 
21), it was found that expressed in degrees, would fall on a 

single curve described as follows: 

Ogjj = 9.195 - 23.734y^ + 60.810yJ - 33.533y 2 3 + 7.391y* - 0. 581y 3 , 


if y £ < 2.5 . 
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( 68 ) 


<* BD = 38.0, if y ^ > 2.5 

where y^ is the distance from apex to centroid of the c g 
distribution from inboard to rj of c g ( max ) measured along the leading 
edge and referred to the half span* This is illustrated in Figure 
8a* 

The progression rate of the breakdown point on delta wings 

at a > a was also analyzed based on Wentz f s data. Although there 
BD 

was considerable scatter in the experimental data (Fig. 8b), a 
single curve based on a least-square analysis was obtained to be 

Ax = 0. 457(Aa) - l-615(Aa) 2 + 0.0303(Aa) 3 - 0,0027(Aa) 4 

+ 0.00009(Aa) 5 , if Aa < 8.0 deg. (69) 

Ax « 0.5392 + 0.0226(Aa), if Aa > 8.0 deg. 
where Aa = a - a,,- and Ax is the nondimensional x-distance, referred 
to the root chord, from the trailing edge to the breakdown point. 

It is known that the vortex strength after breakdown is reduced 
but not vanished. In the VORSTAB code, the sectional c g c at any 
station where vortex breakdown occurs is multiplied by a factor k to 
represent the residual vortex lift. The factor k is determined, 
again, by analyzing Wentz’s data and is found to be (ref. 22) 
k = 0.131 + 0. 384y , if y^ < 1.49 

k - 0.951 - 0,2Q8y + 0.028y 2 , if 1.49 < y ^ < 3.71 

k = 0.5, if y > 3.72 (70) 

Ar 
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In addition, it was found that upper surface slopes in the 
spanwise direction due to thickness distribution will affect the 
movement of the burst point (ref. 23). To model this effect, it is 
assumed that the local angle of attack is changed by an amount equal 
to the spanwise upper surface angle; i.e., tan -1 ( 3z/9y ) . The new 
local angle of attack is used in Equation (69) to determine the 
vortex burst point iteratively because 3z/3y is, in general, not a 
constant. This effect was found to be quite significant for a 
highly cambered wing, such as the F-106B configuration. 

2.3.5 Fuselage Vortex Lift 

If the aspect ratio of a slender thick wing is reduced, 
eventually it becomes a slender body. If the method of suction 
analogy is applicable to the former, it should also be applicable to 
the latter. Based on this assumption, the following method for 
calculating body vortex lift was developed. The method is based on 
the following assumptions and procedures. 

1. As mentioned previously, the attached flow solution is 
obtained with the axial distribution of G. N. Ward’s 
vortex multiplets. 

2. At any axial station, vortex separation starts at a 8- 
circumference location where is minimum and negative. 
This assumption has been shown to be reasonable (ref. 

24). At low angles of attack, may be positive 
everywhere near the nose. In this case, no vortex 
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3 . 


separation is assumed to occur. In reference 24, the 
axial station (x g ) at which the vortex separation starts 
must be assumed or given by experiments. 

At any axial station, the side force component of the 
negative Cp in the region assumed to have vortex 
separation is integrated to produce a sectional side force 
coefficient (see Fig. 9): 


c ~ / r(x)C sin(0)d0 

s ^ r vx; p 


(71) 


The obtained side force is assumed to be the suction force 
produced by the separated vortex. This suction force is 
presumed to be acting at 0 G (Fig. 10) where 


0=0. - A0 (72) 

c min,p 

A0 -he- (73) 

sf 

Equations (72) and (73) imply that the vortex action point 

is located at a distance from 0 . being proportional to 

min, p 

the suction force. Therefore, 

As - rA9 - h r c . (74) 

sf 

Based on Equation (66a), h becomes 


h = 1 / (4t ) (75) 

In reference 25, f is determined to be l//~2 in subsonic 
flow. Therefore, h = 0.5 is used in the program. 

With 0 c calculated from Equation (72), sectional 
normal force coefficient due to the vortex is given by 


C N, v " c sf cos9 c 


(76) 
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Similar to a wing, the augmented body vortex lift exists 
whenever the planview of a body is not of the delta 
type. Therefore, the concept employed in calculating the 
augmented vortex lift for a wing is also applicable for a 
body. Thus, if C g j^ is the total suction force 
coefficient from the nose portion, the augmented normal 
force coefficient y^ is given by 


C N,VA C C sf,N COS0 A^n 


(77) 


where c is the length over which the nose vortex passes 
and Z n is the nose length. Note that 9 a t * le location 
of augmented vortex action point and is assumed to be 
equal to 9 C at the body shoulder. Some calculated results 
for the body vortex lift in both subsonic and supersonic 
flows can be found in Reference 25. 


2.4 Effect of Boundary Layer Separation 

It is known that the classical lifting line theory cannot be 
used for most lifting surfaces of fighter configurations because of 
limitations of the theory on aspect ratios and sweep angles of 
lifting surfaces. Furthermore, existing methods for calculating 
stall and post-stall characteristics of lifting surfaces are not 
satisfactory, particularly for thin wings used on fighter 
configurations. To improve the prediction methodology, a method 
based on utilizing nonlinear section data was developed (ref. 9). 
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In this method, the nonlinear sectional data was used as near- 
field solutions to be matched iteratively with far-field solutions 
obtained from the lifting surface theory# Based on this concept, it 
is assumed that the effect of flow separation is to reduce the local 
angle of attack by Aa at a spanwise station# Therefore, the 
effective angle of attack at any spanwise section becomes 

a = a - a. - a - Aa (78) 

e n i o 

where a n is the geometric angle normal to the section which may have 
dihedral, a^ is the induced angle of attack, a Q is the angle of zero 
lift, and Aa represents a reduction in a n to be calculated due to 
viscous effects. It follows that 

C «(3-D) ' c >. sl " (c, n ' “i _ °o ' 4c0 <?9) 

a 

Assuming c = 2ir/(l - M 2 ) 1 / 2 , Equation (78) can be solved for a^j 
a 


. -1 , C &(3-D ) , 

a. = a - sm [ — — ] - a - Aa 

in c„ o 


(80) 


a 


Let the 2-D sectional lift coefficient evaluated at a_ - 


a^ be 


C t(2-D) and let 


f = C £(2~D) 
C «.(3-D) 


(81) 


Since C ^( 3 _ D ) is computed with an inv.is.cid theory, its value is 
usually larger than c^ 2 -D) if Aa = 0. Therefore, f is usually less 
than 1.0. In this case, a geometric angle of attack (a*) which 
produces the reduced lift can be found. That is, 
sina' = f sina 

n 

or, 
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( 82 ) 


o' = sin ^(f sina ) 
n 

It follows that Aa in Equation (78) becomes 

Aa = a n - a' (83) 

The solution is obtained iteratively as follows: 

1, Assume Aa = 0. 


2. Find a^ from Equation (80). 

3. Calculate f from Equation (81). 

4. Determine Aa from Equations (82) and (83). 

5. Use Aa to reduce a in the 3-D boundary condition to 
determine 

6. Repeat steps 2 through 5 until the successive total lift 
coefficients differ by less than 0.5%. 

For a configuration of moderate aspect ratio with vortex lift, 
the iterative procedures described above are not started until the 
predicted vortex breakdown occurs. Before vortex breakdown, the 
viscous effect is accounted for only once without iteration. The 
latter procedure is also used in calculating lateral-directional 
aerodynamics. That is, lateral-directional aerodynamic 
characteristics are calculated without iteration only after the 
calculation for longitudinal characteristics has converged. In 
addition, when a wing with a separated flow is subjected to a 
positive sideslip or roll rate, part of the left-wing separated 
region may be reduced or suppresssed because of reduced loading. 
Therefore, the reduced lift due to separation may be recovered. 

This effect of lift recovery will make the resulting roll moment 
less negative and is included in the code. 
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Based on experience, this method is applicable to most of the 
current fighter configurations# However, if the lifting surfaces 
are dominated by vortex flow, such as the F-106B configuration, 
inviscid flow should be assumed# 

2. 5 Calculation of Lateral-Directional Stability Derivatives in 
Attached Flow 

The lateral-directional stability derivatives to be discussed 
in this study are g-, p-, and r-derivatives# 


2.5.1 g-Derivatives (C , C , and C« ) 

y 3 8 *3 

The distribution of lifting pressure coefficients in sideslip 
can be written as 


AC = 2(y + y sing) 

P y x 

Therefore, for small sideslip angles. 


(84) 


, 9 S 
\ " ' 9e 


__2 

Sb 


9Y v 

/ + Y x )cos0 y dy 


(85) 


where 9 is the camber angle. 

To compute the wing contribution to C and C n , the leading- 

y g 3 

edge and side edge suction force components due to sideslip must be 
calculated. The side-edge suction force per unit length is given as 
(ref. 13) 


■S (x) = irpG (x) 


( 86 ) 


where G(x) is the singularity parameter of circulation T t (y) at the 
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tip and is related to F t by 


G(x) - (|) 1/2 li»(l -• E f 7 ) 1/2 1% < 87 > 

b 

y^2 

For the antisymmetrical case, the total side-edge suction force can 
be written as 


S (x) = irpCG ± G ) 2 (88) 

t a a 

where the plus sign is for the right wing and the minus sign for the 
left wing, with subscripts a and a denoting symmetrical and 
antisymmetrical cases , respectively. 

The sectional leading-edge thrust coefficient in combined 
symmetrical and antisymmetrical loading can be written as (ref. 13) 


f A 


2 2 
M„ cos A ji 


(c ± .c v 

a a 
cosA„ 


(89) 


where C„ and C are leading-edge singularity parameters, for 
symmetrical and antisymmetrical loading, respectively. If C a due to 
sideslip is denoted by C a and c r due to sideslip only is Ac , then 

P X 3 

the leading-edge thrust coefficient due to sideslip is given by 

2C C 

o o & 3 p 

Ac 88 /I - M cos A„ — (90) 




Therefore, the wing side force due to leading-edge suction in 
sideslip is 

Ac ~ Ac tanA rt (91) 

y e x s * 
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The total wing side force due to sideslip contains four 


sources: 

1. Contribution from the side-edge suction 

The side force along the side edge is 

S = irp(G ± G ) 2 (92) 

t a a 

Therefore, 

9S 9G 

it - irr> <93) 

2. Contribution from the leading-edge suction 

AC = •§■ / c(y)Ac dy (94) 

7 3 S right 7 3 

wing 

3. Contribution from the incremental pressure force due to 
geometric dihedral 

The incremental lifting pressure force, which acts 

normal to a planform, has components in the side force 

direction* The contribution to C can be calculated by 

y 3 

integrating sectional contributions. 

4. Contribution from the induced drag 

The induced drag is assumed to act in the direction 
of the freestream with sideslip. For small sideslip 
angles, the side force due to the induced drag can be 
computed as 

AC = -C 3 

y D < 
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It follows that 


AC 



( 95 ) 


The fuselage contribution to the side force coefficient can be 
computed as 


C 

y (f> 



r G y(f ) dx 


(96) 


where c T is the fuselage sectional side force coefficient* The 

y<f) 

contribution to C„ is obtained by differentiating Equation (96) 
with 6* 

The yawing moment coefficients due to sideslip are obtained by 
taking moment of all the side force components about the z— axis. 


2. 5. 2 p-Derivatives (C , C , and C ) 
r y n l 

P P P 


The roll damping derivative C is computed by integrating the 

Xr 

P 

antisymmetrical pressure force induced by the roll rate multiplied 

by the spanwise moment arm. The side force and yawing moment due to 

roll rate are contributed from the following sources: 

1. Contribution from the leading-edge suction 

The incremental sectional leading-edge thrust 

coefficient due to roll rate is 

2.C C 



TT 2 Z 

— /l - M cos A 
2 00 


l 


a a 


£ 

cosA^ 


( 97 ) 


30 



Therefore , 


2 . 


3. 


Ac = Ac tanA. 
y x Z 

P P 


( 98 ) 


AC V “ e ' / c(y)Ac dy 
right 


(99) 


wing 

Contribution from the side edge suction 


as 3G 

t — - 2 tt pG (± Q a - P -) 
3p M a 3p 


( 100 ) 


Contribution from the incremental pressure distribution 


The concept described in relation to Cy^ is still 


applicable to C 


C can be obtained by taking moment 
P 


of all the side forces due to p about the z~axis. 


2.5.3 r-Derivatives (C , C , and Cn ) 

y r Il r 

The yawing motion can change the pressure distribution on the 
lifting surfaces. The incremental pressure distribution due to yaw 
rate consists of three components: 

1. Due to yaw rate r, a backwash (ry) is produced (Fig. 11) 
where y is the spanwise station. Therefore, 

AC = -%L Y ( 101 ) 

p r v oo y 

2. The sidewash (rx) effect produces 

AC = - -2L y cos <j> (102) 

pr V x 

r 00 

3. Wing-body interaction is calculated by satisfying the 
boundary conditions. 
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The calculation of CL and CL follows the same 

y r n r 


procedures as for computing C and C n • 

y 3 ‘3 


2. 6 Lateral-Directional Aerodynamics in Vortex Flow 

It is indicated in references 26 and 27 that in a sideslip the 
leading-edge vortex on the windward side is pushed inward and 
downward, thus inducing additional vortex lift, while on the leeward 
side it is pushed outward and upward, thus inducing less vortex 
lift* To account for the increase in vortex lift on the windward 
side, a "displacement-type" vortex lift is introduced. In addition, 
vortex breakdown characteristics are also different from those of 
the symmetrical flow situation. The flow mechanism and 
computational methods for sideslip, yawing, and rolling motions are 
described below* 


2,6.1 Vortex Lift of the Displacement Type in Sideslip 

The leading-edge singularity parameter, defined in Equation 
(25), is proportional to the lifting pressure coefficient. As shown 
in Equation (84), due to sideslip, an extra term, Y x sin3, is added 
to the lifting pressure coefficient in symmetric flow. Since the 
vorticity near the leading-edge must be parallel to the leading 
edge, it follows that 

C, , = G (1 ± tanA sin3) (103) 

(3) a l 
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The leading-edge thrust coefficient in sideslip becomes 

/ i _ m 2 2. .1/2 

2 2 2 ^ M oo COS 

c t(6) = (ir/2)C a (± tan V ing + tan V in 0) c ' o ' sX: 


(104) 


2.6.2 Vortex Lift of the Displacement Type in Yawing 


The induced sidewash due to yawing is 


x - x f _ 
g_ = - 2 r 


(105) 


Therefore, the leading-edge singularity parameter due to yawing 
becomes 


C/ \ " C (1 ± 5 tanA ) 
(r) a r i 


(106) 


2.6.3 Vortex Breakdown in Sideslip 


With a sideslip angle B, the angles of attack for vortex 
breakdown on the windward side («br) an< l the leeward side ( otgL^ are 
calculated as follows: 

1. Calculate y from the symmetrical suction distribution 

X/ 

with the leading-edge sweep angles B for the right 

wing and + 6 for the left wing. 

2. Assume that the maximum vortex strength before breakdown 
is unchanged by sideslip. Because of increase in vortex 
lift on the right wing in a positive sideslip, the maximum 
vortex strength would be reached at a lower a (i.e. , obr) 
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than on the right wing. Since the vortex strength is 
represented by c s which is proportional to sin^a, ag R and 
a BL can be obtained by solving the following equation: 

sin^a^p - sin^Og^ + (± 2 sing + tanA^sin^ B)sin^Og R (107) 

where ,f - ff is for ou^. 

Similar to lift recovery in a viscous separated flow 
mentioned earlier, vortex lift recovery is also 
possible. When a slender wing is at an a greater than 
and at a positive sideslip, vortex breakdown on the 
left wing may be suppressed, thereby recovering the 
symmetrical vortex lift to make the rolling moment less 
negative. 

2.6.4 Vortex Breakdown in Yawing 

The effect of yawing on vortex breakdown is similar to that due 
to sideslip. Since the equivalent sideslip (3 r ) is variable along 
the leading edge, the average of 3 r at a given y station and that at 
the tip is used to determine breakdown a. In the program, a BB for 
the yawing motion is determined by a linear interpolation between 
a BD and ot BR for sideslip. 

2.6.5 Vortex Breakdown in Rolling 

Due to a positive roll, the local angle of attack on the right 
wing is increased. It seems that the vortex breakdown angle would 
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be decreased. However, rolling tends to move the centroid of the 
vortex lift distribution outboard (i.e. , to increase y^) . 

Therefore, the vortex breakdown angle is increased. In the present 
code, the vortex breakdown angle for the rolling motion is assumed 
to be equal to a gD for the symmetrical loading. 

3. NUMERICAL RESULTS AND DISCUSSION 

Earlier applications of the methods described herein can be 
found in reference 22 for simple slender wings and in reference 9 
for two complete configurations of airplanes. For simple slender 
wings, the main objective was to show the effect of vortex 
breakdown. On the other hand, for the airplane configurations it 
was imperative to account for the effect of boundary-layer 
separation. For the latter, nonlinear section data including post- 
stall characteristics are required. 

The performance of some theoretical methods for predicting 

airfoil characteristics up to stall has been reviewed by Blascovich 

in reference 28. Rumsey found another method (i.e. the SAAP) to be 

accurate in predicting c* (ref. 29). These methods are typically 

^max 

good for moderate to thick airfoils with angles of attack up to the 
stall angle. At post-stall angles of attack and for thin airfoils 
typically used on fighter wings, these methods are all 
questionable. On a thin airfoil, the flow separation is usually of 
the long bubble type. The latter is not well predicted by existing 
methods. In view of the fact that all thin airfoils tend to behave 
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in the same way as a flat plate (refs. 30 and 31), for the present 
purpose it was decided to modify Eppler's code (ref. 32) by 
including a modified free-wake theory (ref. 31). Typically, a thin 
airfoil will stall at about 8 to 8.5 degrees of angle of attack. 
After stall, the variation of c^, c^, and with angles of attack 

is assumed to be given by the modified free-wake theory. 

To obtain a better understanding in determining how an airplane 
configuration should be modeled for both longitudinal and lateral 
calculation, several existing fighter configurations were 
investigated. They are 

(a) An F/A— 18 configuration with and without leading edge flap 
deflection. 

(b) A generic fighter model. 

(c) An F-16 configuration. 

(d) An F-16XL configuration. 

(e) The F-5 and modified configurations. 

(f) An X-29A forward swept wing airplane. 

(g) An F-106B configuration. 

In the following, some calculated results will be presented. 

(a) An F/A-18 Configuration 

The F/A-18 configuration is relatively complicated in 
geometry. It has highly cambered strakes so that not only is the 
viscous separation effect important on the wing at high angles of 
attack but also the vortex lift is significant. 



The geometric model of the F/A-18 used in the calculation is 
illustrated in Figure 12. Based on calculations with the modified 
Eppler's code, the sectional characteristics shown in Figure 13 were 
used in the VORSTAB calculation. 

In these calculations, it is assumed that there is no boundary- 
layer separation on the strakes. The strake vortex was observed to 
pass through the lower surface of the vertical tails at angles of 
attack lower than 20 degrees (ref. 33). Thus the characteristic 
length for the augmented-vortex lift (c) to be used is negative. 
However, at angles of attack greater than or equal to 20 degrees, 
the augmented-vortex lift on the vertical tails should be 
positive. Forebody vortex lift is assumed to exist. However, 
interaction of the forebody vortices with the vertical tails in the 
lateral— directional motion is not modeled in the present code. No 
additional discrete vortices are placed at the strake-wing juncture 
in the modeling for this configuration or for the F-16 to be 
discussed later. 

Longitudinal characteristics without the leading-edge flap 
deflection are presented in Figure 14. It is seen that is 
slightly underpredicted and C m is more negative for angles of attack 
in the range of 15 to 30 degrees. This is perhaps caused by the 
predicted vortex breakdown being too early in a and viscous 
separation on the wing. Because of a relatively low sweep angle (26 
degrees), the wing is predicted to have vortex breakdown at an angle 
of attack of 12 degrees, instead of merely being displaced upwards 
due to the influence of strake vortices (ref. 35). 
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Lateral-directional coefficients C„ , C 0 , and C„ of the same 

V H n 3 

configuration are presented in Figure 15. At angles of attack below 

15 degrees, the VORSTAB code was capable of predicting well the 

lateral-directional characteristics. However, discrepancy occurs at 

higher angles, most probably due to the strong strake vortex passing 

through the vicinity of the vertical tails. In a positive sideslip, 

the positive side force on the right vertical tail produced by the 

right strake vortex may be greater than the negative side force on 

the left vertical tail produced by the left vortex, thus generating 

a more negative CL and a more positive C 0 in the test results. 
n B *0 

Another possible reason is the effect of wing wake on the vertical 
tails (ref. 36). At present, these types of interference between 
lifting surfaces still can not be accurately predicted by the 


program. 

To see the effect of leading-edge flaps, nonlinear section data 
are needed. The lift data are presented in Figure 16. Longitudinal 
characteristics are well predicted, as shown in Figure 17, except 
the pitching moment coefficients again are more negative at moderate 
angles of attack. Since the leading-edge flap will delay the wing 
separation effect, C n is now more positive due to the reduced 
effect of wing wake on the vertical tails and is shown to be more 


accurately predicted as shown in Figure 18. 


The discrepancy in 


is mainly due to the predicted vortex-breakdown a of the wing in 


sideslip being inaccurate. They are 24.7 and 25.9 degrees for the 


right and left wings, respectively, compared with possible values of 
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20 and 30 degrees in the data. Of course, the strake vortex 
interaction again plays a role in the C Q variation. 

(b) A Generic Fighter Configuration 

This is an airplane model tested in the 12-foot tunnel at NASA 
Langley Research Center (ref. 37). As shown in Figure 19, the model 
includes a cylindrical body with lifting and control surfaces made 
of flat plates. Nonlinear sectional data from reference 31 were 
used in the calculation (Fig. 20). In the calculation for the 
viscous effect, an overrelaxation factor was used for the wing. 
However, the horizontal tail required an underrelaxation factor to 
achieve convergence. This is most probably due to a coplanar wing- 
tail interaction. In this case, not only should the vortex strips 
on the wing and tail be lined up to avoid unrealistic downwash 
induced on the tail, but also the relaxation factor on the tail must 
be reduced to avoid divergence. Based on the experience with other 
configurations, an underrelaxation factor of 0.5 for both surfaces 
should also work well as an initial trial value. Because of a 


relatively long forebody, body vortex lift is assumed to exist in 
the calculation. This latter effect mainly affects the pitching 
moment calculation in the present code. 

Longitudinal characteristics are presented in Figure 21. It is 
seen that the predicted results agree well with data, except that 
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and C* are not correctly predicted. This is probably due to the 

3 

contribution of asymmetric forebody vortices to C n at high a. At a 

3 

positive sideslip, the asymmetric forebody vortices tend to produce 
a positive side force near the nose and hence a more positive C n ^ 
for the whole configuration. The reason for the large change in 
at a > 25 degrees is not precisely known and could be the result of 
interaction between asymmetric forebody vortices with the vertical 
tail. 


3 


(c) An F-16 Configuration 

As seen in Figure 23, this airplane has a slender leading-edge 
extension (or strake). Also, the sectional data (Fig. 24) are used 
only on the wing. Body vortex lift is included in the 
calculation. Like other airplanes which have been modeled, there is 
no difficulty in modeling the lifting surfaces. However, the F-16 
has a distinct inlet which makes the cross sections of an equivalent 
body of revolution change too rapidly. Note that only a body of 
revolution is allowed in the code. In addition, the inlet makes it 
more difficult for the predicted fuselage aerodynamic 
characteristics to be reasonable if the real body shape with the 
nacelle is used. Therefore, smooth body cross sections are 
assumed. As shown in Figure 25, the longitudinal characteristics 
are reasonably well predicted, up to an a of 25°. - , 

Unlike the F/A-18 configuration, the strake for the F-16 model 

was assumed noncambered and untwisted. Therefore, the 

predicted a was lower than that for the F/A-18 (31 vs 33 
BD 
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degrees). This might have affected the lift prediction at high a. 
Another reason for the drop in is that the k-f actor for the 
strake could be too low (see Eq. 70). Note that an inviscid theory 
would overpredict the lift. 

The strake of this configuration is not as large as the one on 
the F/A-18. Therefore, the influence of strake vortices on the 
vertical tail is not as significant. The predicted lateral 
characteristics agree well with the data, as shown in Figure 26. 

(d) An F-16XL Configuration 

This configuration (Fig. 27) consists of the fuselage and the 
horizontal tail of the F-16A and a highly swept cranked wing. As 
described in modeling the F-16A configuration, smooth body cross 
sections are assumed. The inboard and outboard leading-edge swept 
angles are 70 to 50 degrees, respectively. Because the vortex flow 
is expected to be strong, nonlinear sectional data were not used in 
the model. Body vortex lift is not modeled for this configuration. 

Predicted results for longitudinal aerodynamics without viscous 
effect are shown in Figure 28. Again, the prediction agrees well 
with data. When sectional data were specified over the outboard 
portion of the cranked wing, the predicted lift coefficient was too 
low by about 14% at high a (not shown in the figure). However, all 
other characteristics, including the lateral-directional parameters, 
were well predicted. Sideslip derivatives are presented in Figure 
29. Except at a > 30 degrees, both C and C are accurately 

n e 3 

computed. At high a, the 
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discrepancy is most likely caused by the burst vortex flow on the 

right side inducing positive side force on the vertical tail* The 

latter would produce more negative G and positive C Q • The 

n 0 0 

predicted a for vortex breakdown for the inboard portion of the wing 
in a positive sideslip (oijjr) about 32 degrees* 

(e) The F-5 and Modified Configurations 

There are three wing conf igurations in this series of models: 
the basic F-5, an F-5 fuselage with a swept wing, and an F-5 
fuselage with a delta wing (Figs* 30 and 31). The F-5 configuration 
is well known for its long forebody, which can generate strong 
forebody vortices. For strong wing-dominated vortex flow, the 
forebody vortex effect on longitudinal characteristics is minor 
compared with the wing vortex lift. On the other hand, for a 
moderate aspect-ratio wing, the body vortex may induce strong enough 
spanwise pressure gradient on the wing to delay vortex breakdown on 
the latter, in much the same way as a strake does. 

In examining these three configurations, the swept angles of 
the basic and swept wing configurations are not large* Thus they do 
not have strong vortex lift effect and may be subjected to boundary 
separation which is influenced by the forebody vortices. Since 
there is no theoretical method available to determine the degree of 
influence exerted by the forebody vortex on the wing flow field, at 
the present time this can only be done empirically. 

The basic assumption made in this regard is that the forebody 
vortex is to produce a favorable pressure gradient on the wing to 
delay vortex breakdown and the fuselage downstream of the wing 
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leading edge is to cause earlier vortex breakdown. The magnitude of 
these effects is assumed to be proportional to a distance equal to 
the projected area of the corresponding portion of the fuselage on 
the x-y plane divided by the fuselage radius at the strake region 
plus the strake width if it is present. The increase or decrease in 
a B Q is then taken to be the Aa by which a vortex burst point will 
move over that distance. 


The sectional data used in the calculation are again generated 

with the revised Eppler's code and are shown in Figure 32. 

Longitudinal characteristics of the basic configuration are shown in 

Figure 33. The prediction is seen to agree well with the data. 

Again, including the nonlinear section data improves significantly 

the agreement of the lift curve. Sideslip derivatives for both C n ^ 

and C t are also well predicted (Fig. 34). Note that the 

nonlinearity of C with a is mainly predicted by the following 

n g 

mechanisms. A positive sideslip will produce a negative side force 

as usual. However, the freest ream (i.e. V sina) results in a 

negative sideslip to the vertical tail producing a positive side 

force. The wing downwash, viewed as a positive sideslip to the 

vertical tail, will produce a negative side force* In earlier 

versions of the VORSTAB, the wing downwash effect was not included 

so that C n would be decreased rather quickly as a Is increased. 

0 

With both Vasina and wing downwash included, C n ^ will decrease more 
gradually. As a is further increased, the wing vortex will burst 
and the boundary layer will separate to decrease the wing downwash 


at the vertical tail. At the same time, the local dynamic pressure 
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would be lower. All of these effects are to reduce the negative 

side force at the vertical tail. In the present code, an assumption 

is made to remove the wing downwash effect on the vertical tail as 

soon as a is reached. This assumption appears to work well for 
BJJ 

the basic F-5 configuration. At a ) 30 degrees, the forebody 
vortices will have a dominating effect by producing forces on the 
forebody and interference on the vertical tail and are not modeled 
in the code. 

For the F-5 configuration with a swept wing (A = 51 degrees), 
both vortex lift and viscous separation are important. The wing has 
a snag which is known to generate a vortex at high a that produces 
downwash on the outboard wing section and divides the wing vortex 
system into two. Each system is assumed to generate augmented 
vortex lift on the inboard or the outboard section of the wing, as 
the case may be. As shown in Figure 35, the longitudinal 
characteristics are well predicted. The total augmented vortex lift 
is only about 6 percent of the total lift at a = 30 degrees. 

Sideslip derivatives are presented in Figure 36. As shown, 
calculations with or without viscous effect are seen to give similar 
results which agree reasonably well with the data. Note that the 
effect of vortex breakdown is quite important for this configuration 
and was included in both sets of calculations. 

For the configuration with a delta wing, the wing sweep (A = 59 
degrees) is expected to produce a vortex-dominated flow field. 
Therefore, the viscous effects were neglected In the calculation. 

As shown in Figure 37, C L , C D and C m are underestimated 
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at a > 20 degrees (Fig* 37), mainly because a BD i s underpredicted. 

For the sideslip derivatives (Fig. 38), since C is well 

3 

calculated, the discrepancy of must be caused by early unloading 
at the wing tip. It is known that tip loading of a delta wing may 
be reduced by either vortex inboard migration or viscous separation 
because of high loading, or both. 

(f) An X-29A Forward Swept Wing Configuration 

The geometry of the model is presented in Figure 39, The model 
wing is made of a flat plate with a leading-edge radius equal to 
1.3% of the root chord. The control surfaces include the canard, 
wing trailing-edge flaperons, and the fuselage-mounted straked 
flaps. The flaperons and strake flaps were deflected to their full- 
down positions (17.5 degrees and 30 degrees, respectively). The 
calculations are based on the inviscid flow with and without vortex 
breakdown. The deflection angle of the flaperons in the streamwise 
direction was calculated to be 13.6 degrees. It is seen from Figure 
40 that both and Cq are underpredicted perhaps because the 
geometric representation in the calculation is not accurate. All 
necessary geometric quantities for input were measured from Figure 
39. 


In the VORSTAB code, the vortex bursting is implemented to 
progress from the wing tip toward the apex. In the swept-f orward 
wing situation, the leading-edge vortex separation starts at the 
wing tip; and its bursting may start at the wing root and move 
toward the tip, which is opposite to the progression direction 


assumed in the program. In Figure 41, the Cp prediction without 


vortex breakdown appears better than the one with vortex breakdown. 
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(g) An F-106B Configuration 

As shown in Figure 42, the basic F-106B configuration has a 

highly swept wing which is dominated by vortex flow at high angles 

of attack. Therefore, inviscid flow is assumed in the 

calculation. Figure 43 shows good agreement with data for the 

longitudinal characteristics. C D ± s underpredicted at high a 

perhaps because the thrust component of lift acting on the conical 

camber is overpredicted by the code. In the lateral-directional 

calculation (Fig. 44), predicted C„ and C 0 match well with 

n 0 *0 

experimental data up to a = 30 degrees. However, C i s not as well 

y 8 

predicted, probably due to the wing and fuselage interference. 

Figures 45 and 46 show the lateral-directional derivatives due to 

rolling and yawing. The experimental data measured in forced- 

oscillation tests and taken from Reference 32 to compare with the 

predicted C , C , and C« are (C__ + C sina), (C + C since), 
y p n ? *p y p y i n p n g 

and (C^ + C. sina). The test results taken from the same 

P g 

reference to compare with C v , C_ , and C p are (C v - C cosa), 

J r r *r ^r Yg 

~ cosa), and cosa), respectively. This code is 

not designed for unsteady calculation. Therefore, discrepancy at 
high angles of attack is expected. In forced-oscillation tests, 
damping derivatives tend to be higher because of aerodynamic lag 
effect (e.g. viscous and vortex lag). Figure 47 shows the aileron 
control power at an aileron deflection of 7 degrees. The rolling 
moment is mainly contributed by the antisymmetric wing lift due to 
aileron deflection and is correctly computed. The yawing moment is 
not as well predicted at low angles of attack. Figure 48 shows the 
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rudder control power at a rudder deflection of 25 degrees* In this 
case, the yawing moment is more accurately calculated than the 
rolling moment* 

4. CONCLUSIONS 

The VORSTAB code has been shown to be capable of accounting 
for, but not limited to, vortex breakdown effect, viscous effect, 

vortex lift, and suction force, etc. The viscous separation effect 

was included through the use of nonlinear section data* The 

nonlinear section data were calculated with a revised Eppler's code 
for thin airfoils based on a modified free wake theory for post- 
stall characteristics. Calculated results for complete airplane 
configurations showed the longitudinal aerodynamics and the lateral- 
directional derivatives predicted by the VORSTAB code gave 
reasonable results compared with wind-tunnel test data up to an 
angle of attack at which wake interference and forebody vortex 
effect were not significant « Such an angle of attack was found to be 
configuration dependent and might be different for longitudinal and 
lateral-directional aerodynamics. Good results could be expected 
for those conf igurations dominated either by vortex flow or viscous 
separation. Overall, predicted results for longitudinal 
characteristics agreed better with wind-tunnel data than with 
lateral-directional results. One possible reason for the 
disagreement is wake interference between lifting surfaces. Also 
the forebody vortex might transgress the flowfield of the vertical 
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tail leading to an inaccurate prediction of lateral-directional 
derivatives at high angles of attack. In addition, asymmetrical 
forebody vortices in a sideslip would produce significant effect on 
directional derivatives and were not accounted for in the present 
code. 

Based on the experience gained in this investigation, a 
successful modeling for the VORSTAB code would require 

(1) good judgment as to whether a configuration was vortex- 
flow dominated or might require nonlinear section data; 

(2) a good estimate of where the augmented vortex lift might 
occur and the amount of favorable or adverse pressure 
gradient to cause the change in vortex-breakdown angles of 
attack; and 

(3) an accurate geometric representation of a configuration, 
including camber shapes, airfoil geometry, leading-edge 
sweep, etc. 

Typically, forebody and strake vortices were favorable to vortex 
stability of the wing. But the strake vortex would burst earlier 
due to the presence of wing adverse pressure gradient. 

To improve the code, the interference of forebody and strake 
vortices on the vertical tail must be determined. This would 
require the calculation of the locations of these vortices. In 
addition, the effect of wing wake and vortex-burst flow field on the 
empennage would be an important factor in predicting static lateral- 
directional characteristics. For dynamic lateral-directional 
characteristics at high angles of attack, an appropriate unsteady 
aerodynamic theory would be needed. 
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Figure 6. Wing Section for Defining Vortex Action Point (Ref. 11) 



SECTION A- A 

Figure 7. A Control Surface Taken in Cross-Flow Plan (Ref. 11) 
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a Delta Wings 
□ Cropped Delta Wings 
o Cropped Arrow Wings 
— Least Square Average 
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Figure 8a. Vortex Breakdown Angle for Different Planforas (Ref. 22) 




Ao, deg 

Figure 8b. Progression Rate of Vortex Breakdown Points for 
Different Planforms (Ref. 22) 
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Figure 19. A Generic Fighter Model Tested in the 12- ft Tunnel 
at NASA Langley Research Center (Ref. 37) 
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Flat Plate Sectional Data (Ref. 31) 



Figure 20. Sectional Lift Coefficients Used in the Generic 
Fighter Calculation 
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Experimental Data (Ref. 37) 
VORSTAB Calculation 



ANGLE OF ATTACK 

(a) vs a 

Figure 21. Longitudinal Aerodynamic Characteristics for a 
Generic Configuration 
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- 64A204 Airfoil 

- Flat Plate (Ref. 31) 
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ANGLE OF ATTACK 

Figure 24. Sectional Lift Curve for an NACA 64A204 Airfoil 
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Experimental Data (Ref, 38) 

— — VORSTAB Calculation with Viscous Effect 
~ ~ * VORSTAB Calculation without Viscous Effect 


No Experimental Data Available 


9 °- 0.000 


ANGLE OF ATTACK 

Figure 26. Lateral-Directional Derivatives Calculation of an 
F-16 Configuration Based on the Body Axis at 
8 “ 5 Degrees 
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Figure 30. Three-view Sketch of the F-5 Basic Model (Ref. 40) 
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Experimental Data (Ref. 40) 
VORSTAB , with Viscous Effect 

VORSTAB , without Viscous Effect 



Figure 35. Longitudinal Aerodynamic Characteristics of an 
F-5 Swept Wing Configuration 
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Figure 39. Three-view Sketch of the X-29 Model (Ref. 41) 








Experimental Data (Ref . 41) 

VORSTAB Calculation with Vortex Breakdown Effect 
- ' VORSTAB Calculation without Vortex Breakdown Effect 
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(c) C D vs a 


Figure 40. Continued 






Wing area, 1 5.70 ft 2 
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Experimental Data (Ref. 42) 

VORSTAB Calculation without Viscous Effect 



ANGLE OF ATTACK 

(a) vs a 

Figure 43. Longitudinal Aerodynamic Characteristics of an 
F-106B Basic Configuration 
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Forced Oscillation Data (Ref. 42) 

VORSTAB Calculation without Viscous Effect 



!0 25 31 

ANGLE OF ATTACK 


Figure 45. Lateral-Directional Derivatives Calculation of an 

F-106B Configuration Based on the Stability Axis at 

p = 0.02 
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Forced Oscillation Data (Ref. 42) 

VORSTAB Calculation without Viscous Effect 



Figure 46. Lateral-Directional Derivatives Calculation of an 

F-106B Configuration Based on the Stability Axis at 
r = 0.02 
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Experimental Data (Ref. 42) 

VORSTAB Calculation without Viscous Effect 




Figure 47. Calculation of an F-106B Aileron Control Power, 
= 7 Degrees 
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©— Experimental Data (Ref. 42) 

VORSTAB Calculation without Viscous Effect 




Fi-Siire 48. Calculation of an F— 106B Rudder Control Power* 
= 25 Degrees 
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